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ABSTRACT 

Recent observations indicate that some class sources have orders of magnitude higher 
accretion rates than those of class I. We investigated the conditions for the high 
accretion rates of some class sources by numerical calculations, modelling an external 
trigger. For no external trigger, we find that the maximum value of the accretion rate 
is determined by the ratio a of the gravitational energy to the thermal one within 
a flat inner region of the cloud core. The accretion rate reaches ~ 10 _4 MQyr _1 if 
the cloud core has a > 2. For an external trigger we find that the maximum value 
of the accretion rate is proportional to the momentum given to the cloud core. The 
accretion rate reaches > lO _4 M0yr _1 with a momentum of ~ O.lMQkms -1 when 
the initial central density of the cloud core is ~ 10 _18 gcm -3 . A comparison between 
recent observational results for prestellar cores and our no triggered collapse model 
indicates that the flat inner regions of typical prestellar cores are not large enough to 
cause accretion rates of ~ 10 _4 MQyr _1 . Our results show that the triggered collapse 
of the cloud core is more preferable for the origin of the high accretion rates of class 
sources than no triggered collapse. 
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1 INTRODUCTION 

Pre-main-sequence stars are classified according to their 
spectral features. Lada & Wilking (1984) and Lada (1987) 
proposed a classification for pre-main-sequence stars us- 
ing the slope of their infrared Spectral Energy Distribution 
(SED). They insisted that this classification corresponds to 
an evolutional sequence. We state class I, II, and III in 
the order of evolution. In addition to these classifications, 
Andre, Ward-Thompson, & Barsony (1993) proposed a new 
category called class 0. The sub-mm (A > 350/^m) lumi- 
nosity of class sources is 5 x 10 -3 times higher than to- 
tal luminosities, which indicates that class sources are 
surrounded by significantly large amounts of circumstellar 
material. Andre & Montmerle (1994) insisted that class 
sources are younger than those of class I. 

The accretion rates of class sources have been believed 
to be higher than those of class I by comparing the obser- 
vations and theories. Bontemps et al. (1996) observed the 
CO outflow activities of YSOs nearby regions of star forma- 
tion, and revealed that some class sources have an order 
of magnitude larger momentum flux than those of class I. 
Most ejection models predict that the momentum flux of 
the outflow is proportional to the accretion rate (Pelletier 
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& Pudritz 1992; Wardle & Konigl 1993; Shu et al. 1994). 
This means that these class sources have higher accre- 
tion rates than those of class I. Moreover, Jayawardhana, 
Hartmann, & Calvet (2001) carried out radiative transfer 
calculations of infalling, dusty envelopes surrounding em- 
bedded protostars. They concluded that the prototype class 
source VLA 1623 in the p Oph cloud must have an ac- 
cretion rate > 10 _4 MQyr _1 by matching the far-infrared 
peak in the SED. This value is much higher than the typical 
accretion rates of ~ 4 x 10 _6 MQyr _1 of the class I sources 
in Taurus- Auriga, whose values were obtained by similar ra- 
diative transfer calculations (Kenyon, Calvet, & Hartmann 
1993). Recently, Di Francesco et al. (2001) observed optically 
thick line H 2 CO(3i 2 - 2u) in the class source NGC1333 
IRAS 4A, and detected the inverse P-Cygni profile in the 
line. They derived a high infall velocity and a high accretion 
rate of 1.1 x lO" 4 M yr _1 . 

One possibility is that the high accretion rates of class 
sources are related to the initial density profiles before 
gravitational collapses. If the starting point of star forma- 
tion is a singular isothermal sphere (e.g. Shu, Adams, & 



Lizano 1987), the accretion rate is constant at 



l/G 



through the entire accretion phase, where c s and G de- 
note the sound speed and the gravitational constant, re- 
spectively. The time variation of the accretion rate can be 
caused by an initial density profile different from that of a 
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singular isothermal sphere (Foster & Chevalier 1993; Hen- 
riksen, Andre & Bontemps 1997 hereafter HAB; Whitworth 
& Ward- Thompson 2001). HAB reproduced high accretion 
rates of class sources by assuming that a prestellar core 
has a fiat inner region surrounded by a power-law envelope. 
In their model, the accretion rate is high in the earlier ac- 
cretion phase, and declines later. Their model is consistent 
with the picture that the accretion rates of class sources 
are higher than those of class I. 

In this paper we investigate quantitatively the condi- 
tions for the high accretion rates of class sources based 
on hydrodynamic calculations. When we take into account 
the pressure effect, even though the analytic HAB's model 
neglects it, we must assume that the flat inner regions is 
more massive in order to cause high accretion rates. Recent 
observations reveal both the masses and sizes of the flat in- 
ner regions of prestellar cores (Ward- Thompson, Motte, and 
Andre 1999). Assuming a gas temperature of 10 K, the flat 
inner regions of typical prestellar cores have thermal ener- 
gies that are comparable to the gravitational ones. It seems 
that such cores do not cause high accretion rates. 

We have investigated collapses with and without an ex- 
ternal trigger. If star formation is triggered by an interstellar 
shock wave, a high accretion rate is expected due to com- 
pression by the shock wave. The cloud cores interact with 
the shock wave under various circumstances. For instance, 
triggered star formation by a distant supernova is suggested 
in the Upper Sco association (Walter et al. 1994;Preibisch 
& Zinnecker 1999) and the p Oph cloud (Vrba 1977;Loren 
& Wootten 1986;de Geus 1992;Preibisch & Zinnecker 1999). 
The star formation in NGC 2265 IRS in the Cone Nebula 
seems to be triggered by stellar wind from a B2 star (Thomp- 
son et al. 1998). The triggered formation scenario is also 
proposed for the origin of the solar system based on studies 
of meteorites (Cameron & Truran 1977). The presence of 
short-lived isotopes in meteorites, which should be provided 
by stellar ejecta from a supernova, or an asymptotic giant 
branch star, leads to this scenario. 

Many numerical simulations have been carried out to 
investigate the triggered collapse scenario (Boss 1995; Fos- 
ter & Boss 1996; Vanhala & Cameron 1998). Hennebelle 
et al. (2002) investigated the evolutions of a Bonnor-Ebert 
sphere which is subjected to an increase in the external pres- 
sure. They demonstrated that the features of the density 
and velocity field of prestellar cores and protostars are well 
reproduced, and that the accretion rate is high in the ear- 
lier accretion phase. In this paper, we examine the possibil- 
ity that an external trigger causes a high accretion rate of 
~ 10~ 4 MQyr -1 . We investigate the conditions for the high 
accretion rates of class sources, pursuing numerical calcu- 
lations of the collapse of the cloud core triggered impulsively 
by an external shock wave. 

In Section [5] we describe our models and numerical 
method. In Section |3] we present the numerical results. In 
Section^] we summarize the main results and discuss them. 



2 MODELS AND NUMERICAL METHOD 
2.1 Basic equations and models 

The fluid equations in spherical coordinates are 



dv , 2 dP Gm(r) 
— = -4nr 



dt 



dm 



(1) 



(2) 



where p is the density, v the radial velocity, P the pressure, 
and m the mass within radius r, 



m(r) = / 47T7-' 2 p(r')dr' 



(3) 



Although the density of the cloud core increases drastically 
in the process of star formation, radiative cooling is effec- 
tive except for the high-density region in the centre (Larson 
1969). We assume that the medium is an isothermal gas, 



P = c s p. 



(4) 



The gas temperature T is set at 10 K (c s = 215ms" 1 ), 
which is the typical temperature for globules and cloud cores 
(e.g. Clemens, Yun, and Heyer 1991) in our calculations. 
The mean molecular weight and the ratio of the specific 
heats are set at 2.5 (assuming cosmic abundance) and 7/5, 
respectively. 

We adopt a Plummer-like density profile as an ini- 
tial condition (Whitworth & Ward-Thompson 2001). Recent 
high angular resolution observations of the precollapse cloud 
support this Plummer-like density profile. Ward- Thompson 
et al. (1994) demonstrated that prestellar cores have flat in- 
ner regions surrounded by power-law envelopes based on the 
observation of starless NH3 cores. We assume the following 
initial density profile: 



P0(r) = Pflat 



Rfi. 



(5) 



At a radius smaller than Rfi a t, the density is nearly con- 
stant, pfiat- At a radius larger than Rfi a t, the density is 
proportional to r~ v . We fix the parameter r\ = 2 in our 
models presented in section 3. 



2.2 Numerical methods 

We solved numerically equations £0-01 using a finite differ- 
ence scheme. These equations were integrated by Godunov's 
method. This Lagrangian hydrodynamical code has a second 
order accuracy in both space and time. We set the inner 
(r = 0) boundary conditions such that P and p have a zero 
gradient, and v = 0, and the outer boundary conditions that 
P, p, and v have a zero gradient. We used the non-uniform 
grids, Arrii+i = (1 + t)Am„ where Arrii is the mass inside 
the ith grid from the centre. The value of e was set at 3 x 10 -3 
in order to obtain a high resolution at the centre. We used 
6000 grids in our calculations. For the cloud core itself of 
the mass 3 Mq, 2737 grids were allocated. Outside of the 
cloud core, the density is equal to that of core edge, and the 
gravitational force was excluded from our calculation. 

We used the sink-cell method adopted by Boss & Black 
(1982) to calculate the later accretion phase. In many nu- 
merical simulations which treated the collapse of the isother- 
mal sphere, the density and the infall velocity of the central 
region increase drastically during a runaway collapse phase. 
It is difficult to continue calculations to a later stage, because 
the time step that satisfies the CFL condition is too small. 
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Table 1. Model parameters and the maximum value of accretion 
rate M max for the no triggered collapse model. 



0.5 

time (10 s yr) 

Figure 1. Accretion rate as a function of time. The solid line 
denotes the analytic solution given by Whitworth & Ward- 
Thompson (2001). The circles and the triangles denote the case 
of T =10 K and T = 0.5 K, respectively. 



When the central density reaches a reference density value, 
psink, cells within radius r ain ^ are treated as sink-cells: the 
accreted gas into sink-cells is treated as point like mass at 
the centre. Since the infall velocity near the sink-cell exceeds 
the sound velocity, the effect of the sink-cell does not travel 
to the outer gas. 

As a test for the numerical code, we compared our result 
with the analytic solution of Whitworth & Ward- Thompson 
(2001). We evaluated the accretion rate at the radius r acc = 
300 AU. Fig. Q shows comparison between the analytic so- 
lution and our results. Each parameter was set as follows: 
Pfiat = 3.0 x TO"" 18 gem' 3 , Rfi at = 5350 AU, and r) = 4. 
The gas temperatures are set at 10 K and 0.5 K. The re- 
sult in the case of T = 0.5 K agrees well with the analytic 
solution. On the other hand, in the case of T = 10 K, the 
maximum value of accretion rate is less than that of the 
analytic solution because of a pressure effect. 



3 RESULTS 

We studied two different models in order to seek for the 
conditions for a high accretion rate of ~ 10 -4 MQyr -1 . In 
Section ^. H and Section[<0] we consider the collapse of cloud 
core without and with an external trigger, respectively. We 
call the former the no triggered collapse model, and the lat- 
ter triggered collapse model in this paper. 



3.1 No triggered collapse model 

Here, we present the numerical results of the no triggered 
collapse model. We summarize the model parameters and 
the maximum value of accretion rate in Tabled Fig. [5] and 
Fig. |3] show the evolutions of the density and the velocity 
profiles in model Al. Here, we use the time r = t/tff, 
which is normalized by the initial free fall time, tff = 

I 32 Gp fl t j > a t t ne centre. In the runaway collapse phase 
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Figure 2. Evolution of the density profile in model Al. The 
number attached to each line denotes the time normalized by 
the initial free fall time at the centre. 
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Figure 3. Evolution of the velocity profile in model Al. The 
number attached to each line is the same as in Fig. |5] 



4 Kazutaka Motoyama and Tatsuo Yoshida 



(t < 3.09), the flat inner region shrinks and its density in- 
creases as the collapse proceeds. The density is proportional 
to r~ 2 in the outer envelope. The infall velocity is nearly 
uniform at ~ 3.3 c s in the outer region. After core formation 
(r > 3.09), the density and the infall velocity are propor- 
tional to r _3//2 and r _1,/2 in the central region, respectively. 

Fig.^Jshows the accretion rates of each model as a func- 
tion of the envelope mass, which is used as an evolutionary 
indicator. The envelope mass denotes the mass that remains 
outside of radius r acc . In all models, the accretion rates in- 
crease rapidly during the earlier collapse phases. After the 
accretion rate reaches the maximum value, M max , it gradu- 
ally declines. The maximum value of accretion rate increases 
as the ratio a of the gravitational energy to the thermal one 
within the inner flat region is large. Fig. ^indicates that the 
value of a is larger than 2 in order to cause an accretion 
rate of ~ 10 _4 MQyr _1 . We compared our results with the 
accretion rates of the pressure-free models. For models with 
a > 2, the maximum values of the accretion rates can be 
fitted well with pressure-free solutions. 

We can estimate the accretion rate for a larger a, ne- 
glecting pressure effect. The inner flat region with mass 
Mfiat collapses within 



M = t ff (R flat )-t 



ff 



V3(4 - V) 



tff(Rflat), (6) 



where tff(Rfi a t) represents the free fall time of the mean 
density, p(Rfiat) = [4-Kpfi a tR 3 f i at (l ~ f )]/(§ irR) lat ), within 
the radius Rfi at . The accretion rate is given by 
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Figure 4. Accretion rate as a function of the envelope mass. 
The open circles, open triangles, open squares, and filled squares 
denote model Al, Bl, Dl, and D2, respectively. The number at- 
tached to each plot denotes the ratio a of the gravitational energy 
to the thermal one within the inner flat region. 



Table 2. Model parameters, the arrival times at the centre of 
shock waves t s , and the maximum values of the accretion rate 
for triggered collapse models. The value of the initial free fall 
time at the centre tff is 6.65 X 10 4 yr in the case of pfi a ± = 



1.0 X 10~ iB gcm- 



M 



Mflai 

At 



7T (l 



2 



( 2GM 3 flat 



1/2 



(T) 
(8) 



If we consider the inner flat region as approximately an uni- 
form sphere, the ratio a is estimated to be 

a ~ 2GM flat /(5c 2 s R flat ) 



0.42 



Mfiat \ ( Rfiat y 1 ( t y 1 

O.ISMqJ V3400AU/ U0K7 ' { > 



3400 AU, 

Substituting this relation into equation ® gives 
T 



A I 



4.24 x 10" 5 a 3 / 2 (y^) 3/2 MQyr" 1 



(10) 



This estimation neglecting pressure effect is appropriate for 
the cases with a > 2. Then, the accretion rates are essen- 
tially determined by Rfi at and Mfi at , although equation 
HOII seems to depend on T. A similar estimation is also 
found in Ogino, Tomisaka, & Nakamura (1999). 



3.2 Triggered collapse model 

Here, we present numerical results of the triggered collapse 
model. We give some grids outside of the cloud core the 
velocity to a centre, Vi„i . The total amount of mass in these 
grids is Mini. A similar method is used in Boss (1995). We 
summarize the model parameters, the arrival times at the 
centre of shock waves t s , and the maximum values of the 
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accretion rates in Tabled We adopt the density distribution 
of model Al as an initial condition, of which a is the smallest 
in no triggered collapse models. 

Fig.|H|and Fig.|S|show the evolutions of the density and 
the velocity profiles in model E20, respectively. While the 
shock front propagates to the centre (r < 2.14), the shock 
wave compresses the envelope gas. As a result, the cloud core 
has a centrally concentrated configuration relative to the ini- 
tial state. The density of this region is approximately uni- 
form and one order higher than the density of the preshock 
region, where the shock wave has not yet arrived, and the 
density profile of the outer region is steeper than oc r~ 2 . 
The propagating speed of the shock front is nearly constant 
during compression. The velocity of the preshock region is 
nearly at rest, and the velocity of the postshock region is 
uniform. Compared with no triggered model, the envelopes 
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are denser and have steeper outer-regions during the earlier 
collapse phases, which correspond to the class phases. 

Fig. represents the accretion rate in model E20 com- 
pared with model Al. In model E20, the accretion rate takes 
the maximum value (r = 2.24) just after the collapse of 
the flat inner region. After that, the accretion rate declines 
rapidly. In the later phase there is no large difference be- 
tween the accretion rates of model Al and model E20. Other 
models have the same time evolution as model E20. 

Our results suggest that the maximum values of the 
accretion rates are proportional to momenta given to the 
cloud cores. Fig.|H]plots the maximum values of the accretion 
rates versus the momenta given to the cloud cores, Mi n iVi n i- 
The models with the same momenta show similar values of 
M mol and t s : i.e., E10 and F5, E20 and F10. The maximum 
value of accretion rate is determined by the momentum given 
to the cloud core. A momentum > 0.1MQkms~ 1 causes an 
accretion rate of > 10 MqJT . 

The maximum value of the accretion rate can be es- 
timated as follows. From equation |), a of the centrally 
concentrated flat region compressed by external shock wave 
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Figure 5. Evolution of the density profile in model E20. The 
number attached to each line is the same as in Fig. |2] 



2G M cor 

<jl, s ±lcor, 



(11) 



Here, M com and R CO m denote the mass and the radius of 
this compressed region, respectively. If the arrival time at 
the centre of the shock wave t 3 is shorter than the time tnt, 
when the central density reaches p S i n k, for the no triggered 
case (i.e. model Al), we can assume that the density and 
the velocity of the preshock region remain unchanged from 
the initial state. Using the relation of the compression ratio 
for isothermal shock, we obtain p com /pfi a t ~ {vcom/c s ) 2 , 
where p cora and v com denote the density and the velocity of 
the compressed region, respectively. Thus, equation Hill is 
rewritten as 



(12) 



2G ( hrpfiat y/ 3 
° ~ 5ci V 3ci J 

The momentum within the compressed flat region is con- 
served if the impulse I pre given from the preshock region is 
negligible. This impulse 

Ipre = / 4-Kr 2 c 2 s p (r s )dt (13) 
Jo 

is negligible compared with Mi n iVi„i, assuming a constant 
velocity of the shock front (i.e. r s = -R c (l — t/t s )), where r s 
and R c — 1.7 x 10 4 AU denotes the distance from the centre 
to the shock front and the radius of the cloud core, respec- 
tively. We can replace M com v com with Mi„iVi n i in equation 
11211 . Substituting this relation into equation HIGH gives 
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As shown Fig. |5] our results are classified into three cat- 
egories according to the comparison between the arrival 
time at the centre of the shock wave t s and the required 
time for the central density to reach p 3 i n k, tnt = 3.11 1//. 
In the case of t s -C t n t (Mi ni v in i > 0.3MQkms _1 ), our 
numerical results agrees well with the relation of equation 
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Figure 6. Evolution of the velocity profile in model E20. The 
number attached to each line is the same as in Fig. H\ 



1141 . In the case of t s ~ t n t (0.05MQkms _1 < Mi„iVi„i < 
0.3MQkms _1 ), the maximum values of the accretion rates 
are higher than the relation of equation 114H . because there 
are effects of both compression by the shock wave and con- 
traction by self-gravity. In the case of t s > tnt {MiniVini < 
0.05MQkms _1 ), the effect of the compression by the shock 
wave is small. 



4 SUMMARY AND DISCUSSION 

We summarize our main results as follows: 

(i) The collapse of cloud core with a > 2 can be treated as 
pressure-free, and the maximum value of the accretion rate 
is proportional to a 3 / 2 in the no triggered collapse, where a 
is the ratio of the gravitational energy to the thermal one 
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Figure 7. Accretion rates as a function of the envelope mass. 
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E20, respectively. 




0.2 0.4 0.6 

momentum (M© km s ) 

Figure 8. Maximum values of the accretion rates M ma x versus 
the momenta given to cloud cores, Mi n iVi n i. The double open 
circle denotes model Al. The open diamond denotes model D3. 
The open triangles denote models E5, E10, E15, and E20. The 
open circles denote models F5, F10, F15, and F20. The open 
squares denote models G10 and G15. The solid line denotes the 
relation of equation 1141 . 



within the flat inner region of the cloud core. A flat inner re- 
gion with a > 2 causes an accretion rate of > lCP 4 MQyr -1 . 

(ii) For triggered collapse, the cloud core has a centrally 
concentrated configuration relative to the initial state as the 
result of compression by the shock wave. Just after the shock 
front reaches the centre, when the cloud core has compressed 
the flat inner region and outer steep envelope, the accretion 
rate attains the maximum value. 

(iii) The maximum value of the accretion rate is propor- 
tional to the momentum given to the cloud core in the trig- 
gered collapse. The momentum > O.lMQkms -1 causes an 



Hennebelle et al. (2003) investigated the evolutions of 
a Bonnor-Ebert sphere which is subjected to an increase in 
the external pressure with Smoothed Particle Hydrodynam- 
ics simulations. Their simulations indicate that the prestellar 
core has a flat inner region and a steep outer envelope, and 
approximately a uniform infall velocity profile. Their results 
for the fast compression cases, cf> < 1, are similar to ours 
concerning the density and the velocity profiles, where their 
parameter (j> is defined as the ratio of the time scale on which 
the external pressure doubles to the initial sound cross- 
ing time. These features are consistent with observations of 
prestellar cores. However, in their calculations the maximum 
value of the accretion rate reaches only ~ 3 x 10 _5 MQyr _1 . 
This low accretion rate is due to their low initial central 
density, pfiat ~ 6 x 10 -20 gcm -3 , which is two orders lower 



than that of our model, pfiat = 10 18 gcm 3 . 

We discuss the masses and the sizes of the flat inner 
regions of prestellar cores. Although in the p Oph cloud 
there is a high-accretion class source, i.e. VLA 1623, a 
comparison between the observational results and our no 
triggered collapse model indicates that the flat inner regions 
of prestellar cores are not large enough to cause accretion 
rates of > lO _4 M0yr _1 . The mass and the radius of the flat 
inner region of the prestellar core L1689B in the p Oph com- 
plex are estimated to be ~ 0.33 Mq and ~ 4000AU (Andre, 
Ward- Thompson, & Motte 1996). The parameters of this 
prestellar core nearly coincide with those of model Bl, in 
which the accretion rate and a are 5.45 x 10 MQyr" 1 and 
0.67, respectively. Moreover, Ward- Thompson et al. (1999) 
recently observed eight isolated prestellar cores. They es- 
timated the typical mass and the radius of the flat inner 
regions to be ~ 0.7 Mq and ~ 4000AU. The parameters of 
this prestellar core nearly coincide with those of model C, 
in which the accretion rate and a are 9.77 x lO -5 M0yr -1 
and 1.3, respectively. These results indicate the difficulty of 
evolution from the typical prestellar core into the observed 
high- accretion class source without an external trigger. 

In regions where stars form in cluster, some observa- 
tional results indicate the triggered collapse scenario. Tachi- 
hara et al. (2002) found that clouds with cluster formation 
have head-tail structures, studying statistically 179 cloud 
cores in nearby star-forming regions. Such structures indi- 
cate interactions with the external shock wave. Moreover, 
the class envelopes in Perseus are 3 to 12 times denser 
than singular isothermal sphere at T = 10 K (Motte and 
Andre 2001). The high density of class sources can be in- 
terpreted as being the results of compression by an external 
shock wave. These facts suggest that star formations are 
caused by external triggers in these regions. 

Momentum plays an important role in triggered star 
formation. Foster & Boss (1996) investigated the interac- 
tion of a plane wave and the cloud core by two-dimensional 
calculations. They showed that a fast adiabatic shock wave 
is liable to disrupt the cloud core, and that a slow isother- 
mal shock wave causes a collapse of the cloud core. They 
also showed that the momentum given to the cloud core is 
more essential than the kinetic energy in predicting whether 
induced collapse occurs or not. 
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Next, we discuss the magnitude of the momentum in 
triggered collapse. Foster & Boss (1996) estimated the mo- 
mentum given to the cloud core by an explosion with an 
initial kinetic energy of E = 2 x 10 50 erg. They assumed 
that the evolution of a remnant switches from the Sedov- 
Taylor expansion phase to an isothermal expansion phase at 
a shell temperature of 10 6 K. They adopted a mean molec- 
ular weight of fi = 0.61. This switch occurs at a radius 
of ~ 14pc(E/2 x 10 50 erg) 1/3 (n/lcm" 3 )" 1/3 . The momen- 
tum impacting on the cloud core at a distance R = 29 pc 
from a supernova is estimated to be ~ 0.1 MQkms -1 (E/2 x 
10 50 erg)(R c /1.7 x 10 4 AU) 2 (R/29 pc)~ 2 , where R c is the ra- 
dius of the cloud core in our triggered collapse model. There- 
fore, a supernova located up to 29 pc away from cloud cores 
would cause a high accretion rate of > 10~ 4 MQyr" 1 . Al- 
though the cloud core is impacted by a planar shock in 
this situation, our calculation is spherically symmetric. If 
the pressure of the post-shock region is sufficiently compa- 
rable to that of the cloud core, and the shock wave does 
not significantly compress before he shock front passes the 
core (i.e. t s > 2R c /vi ni ), the cloud core would be spherically 
compressed. Otherwise, the cloud core would be compressed 
from one side (Lim, Hartquist, & Williams 2001). The ac- 
cretion rate would have a somewhat lower value than in the 
case of spherically symmetry. In the p Oph cloud scenarios 
of shock triggered star formation suggested by many authors 
(Vrba 1977; Loren & Wootten 1986; de Geus 1992; Preibisch 
& Zinnecker 1999), an expanding shell of Hi gas created by 
a supernova in the Sco OB2 association induces star forma- 
tion in the p Oph cloud at a distance of about 20 pc from 
the Upper-Scorpius subgroup of the Sco OB2. 

We have investigated the conditions for the high accre- 
tion rates of class sources. Our results indicate that it is 
difficult to cause a high accretion rate of > 10 _4 MQyr^ 1 
with no triggered collapses of typical prestellar cores. For 
no triggered collapse, the flat inner region of prestellar core 
must be a > 2 in order to cause this high accretion rate. 
Otherwise, an external trigger is necessary for a high accre- 
tion rate. Future high-resolution observations (e.g. ALMA) 
will reveal the detailed density and velocity structures of 
prestellar and protostellar cores. Observations with high ve- 
locity resolution give accurate accretion rates of protostars 
from their line profiles. If we know the detailed structures 
of the envelopes and accretion rates of class sources simul- 
taneously, it would reveal whether an external trigger affect 
the accretion rate. 
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